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Abstract 

We show that the most general two-matrix model with bilinear coupling under- 
lies c = 1 string theory. More precisely we prove that W^i+oo constraints, a subset 
of the correlation functions and the integrable hierarchy characterizing such two- 
matrix model, correspond exactly to the W\ +OQ constraints, to the discrete tachyon 
correlation functions and to the integrable hierarchy of the c = 1 string theory. 



1 Introduction 

c = 1 string theory is a well-studied subject. The spectrum of the states is well- 
known, the correlation functions of the so-called discrete tachyons have been calcu- 
lated at least in low genus, the underlying W constraints have been written down in 
explicit form (for reviews on the subject, see for example M). Recently the attention 
of several people has shifted to the problem of identifying the integrable hierarchy 
underlying c = 1 string theory and the topological field theory it may be related 
to. The hierarchy proposed f2|,|| is the dispersionless Toda lattice hierarchy. A 
Landau-Ginzburg potential has also been identified 0],|5[,[@] > @ for this theory. 

The purpose of this letter is to show that the two-matrix model with bilinear 
coupling H provide a unified framework for all these scattered elements of c = 1 the 
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theory. In fact in section 2, after recalling the main features of two-matrix models, 
we show that the W^i+oo constraints proposed in order to interpret the c = 1 string 
theory are nothing but the V^i+oo constraints of two-matrix models. Using the latter 
in section 3 we show how to calculate correlation functions of discrete tachyons and 
verify that they coincide at least in genus with the ones in the literature. Section 
4 is devoted to the comparison between the hierarchies: it was found in || that 
two-matrix models are characterized by the Toda lattice hierarchy; we show that in 
genus this is exactly the integrable structure suggested in Q and in [|| . We show in 
particular that some ad hoc assumptions made in the latter references are a natural 
consequence of the two-matrix model structure. 

2 The two— matrix model and its VFi +00 con- 
straints 

Two-matrix models are defined by the partition function 

Z N (t,c) = J dM 1 dM 2 e TrU 

where M\ and M 2 are Hermitian N x N matrices and 

U = V 1 + V 2 +gM 1 M 2 

with potentials 

Pa 

V a = Y,t^rM r a a = 1,2. (2.1) 
r=l 

The paS are positive integers, characterizing the particular model. In this paper we 
consider p\ and p 2 as arbitrarily large numbers and denote them simply as oo. In 
the following we refer to this model as the two-matrix model. 

We showed in || that one can map this functional integral problem into a linear 
integrable system together with definite coupling constraints. For the sake of clarity 
and with the purpose of introducing formulas and notation we will need later, we 
review this part once again. The ordinary procedure to calculate the partition func- 
tion consists of three steps P| , PHI : (i) one integrates out the angular parts so that 
only the integrations over the eigenvalues are left; (ii) one introduces the orthogonal 
polynomials 

£n.(Ai) = A™ + lower powers, Vn(^2) = A 2 + lower powers 

which satisfy the orthogonality relations 

dA 1 dA 2 e„(A 1 )e yi ( Al ) +v ' 2 ( A2 ) + f AlA2 7 ?m (A 2 ) = h n (t,g)5 nm (2.2) 



(Hi), using the orthogonality relation (|2.2j ) and the properties of the Vandermonde 
determinants, one can easily calculate the partition function 

N-l 

Z N (t,g) = const N\ ]J hi (2.3) 

i=0 
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Knowing the partition function means knowing the coefficients h n (t,gys. 

The information concerning the latter can be encoded in a suitable linear system 
plus some coupling conditions, together with the reconstruction formulas for Zj<j. 
But before we pass to that we need some convenient notations. For any matrix M, 
we define 

{M)ij= M i:j ^-, Mij = M jU M t (j) = Mjj-i. 

As usual we introduce the natural gradation 

deg[Eij] =j-i 

and, for any given matrix M, if all its non-zero elements have degrees in the interval 
[a, b], then we will simply write: M E [a,b]. Moreover M+ will denote the upper 
triangular part of M (including the main diagonal), while M_ = M — M + . We will 
write 

N-x 

Tr(M) = Y, M a 

i=0 

Let us come now to the above mentioned linear system and coupling conditions. 
First it is convenient to pass from the basis of orthogonal polynomials to the basis 
of orthogonal functions 

*n(Ai) = e v ^kn(Xi), $n(A 2 ) = e v ^ Vn (\ 2 ). 
The orthogonality relation ( |2.2| ) becomes 

J d\ 1 d\ 2 y n (\i)e 9XlX2 <f>m(\2) = S nm h n (t,g). (2.4) 

As usual we will denote the semi-infinite column vectors with components ^o, ^2, • • • , 
and $0; 3*1 > ^2; • • • , by \P and <&, respectively. 

Next we introduce the following Q-type matrices 

J dX 1 dX 2 ^! n {X l )X cl e^ x ^ m {X 2 ) = Q nm (a)h m = Q mn (a)h n , a = 1,2. (2.5) 

Both Q(l) and Q(2) are Jacobi matrices: their pure upper triangular part is J+ = 

Beside the above Q matrices, we will need two P-type matrices, defined by 

J dX^X 2 (^-tfn(Ai)) e 9AlA2 $ m (A 2 ) = P nm {l)h m (2.6) 

J dX 1 d\ 2 ^ n {X 1 )e 9XlX ^—^ m (X 2 ))=P mn {2)h n (2.7) 

The following relations hold 

P(1)+ 5 Q(2) = 0, gQ(l)+P(2) = 0, (2.8) 
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It is just these coupling conditions that lead to the famous Wi+oo-constraints on the 
partition function. From them it follows at once that 

Q(l)G[-oo,l], Q(2) e [—1, oo] 

The derivation of the linear systems associated to the two-matrix model is very 
simple. We take the derivatives of eqs.( |2.4D with respect to the time parameters t ar , 
and use eqs.( |2.5| ). We get in this way the time evolution of \I/ and the first discrete 
linear system: 



f Q(l)tf(Ai) = Ai*(Ai), 



3<M>(Ai) = Q% (l)*(Ai), 



c)t. 



-*(Ai 



(2.9) 



^*(Ai) = P(l)*(A x ). 
The corresponding consistency conditions are 

1 



[0(1), ^(i)] 
^-0(1) = [Q(i), &(<*)], 

<} p(i) = [g*(i),p(i)], 
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Of 



(2.10a) 
(2.10b) 

P(1) = [P(1), Q k _{2)] (2.10c) 



where a = 1, 2. 

By studying the evolution of $ we get the second discrete linear system 
Q(2)*(A 2 ) = A 2 $(A 2 ), 
$(A 2 ) = Q£(2)*(A 2 ), 



2.*- 

^$(A 2 ) = -Q fc _(l)^(A 2 ), 

[ &*(A 2 ) = P(2)¥(A 2 ). 
The corresponding consistency conditions are 

[Q(2), P(2)] = l 

Q( 2 ) = [Q(2), Q fc _(a)], 



(2.11) 



P(2) = [P(2), Q fc _(l)], 







dt 



2.k 



(2.12a) 
(2.12b) 

P(2) = [Q k _(2),P{2)\ (2.12c) 



where a = 1,2. Eqs. (|2.10b|) , (|2.10d ),( ^T2b|) , (|2.12cD define the Toda lattice hierarchy. 

The third element we need is the link between the quantities that appear in the 
linear systems and in the coupling conditions with the original partition function. 
We have 



-lnZ N (t,g)=Ti(Q r (a)), a = 1,2 



(2.13) 
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It is evident that, by using the flow equations above we can express all the derivatives 
of Ztv in terms of the elements of the Q matrices. For example 



-lnZ N (t,g) = (Q r (a)) NN _ i , a = 1,2 (2.14) 



dt^idtc 

Knowing all the derivatives with respect to the coupling parameters we can recon- 
struct the partition function up to an overall integration constant. 

We also remark that, since the RHS of the above equations is always defined, 
they give us a definition of Zn even in subsets of the parameter space where the 
path-integral is ill-defined. 

Another consequence of eq.( |2.13| ) and of the definitions of Q{a) and ^, is that we 
can write 

oo °° 1 (9 

ln^iv = £ t r X[ + MnAi - £ — — lnZ^ (2.15) 

r=l r=l 1 ' 

from which we have 

P(l) = £ rtrQilJ- 1 + NQ{1)- 1 + £ Q{iyr-l2^E (2.16) 

r=l r=l r 

Similar formulas hold for ln$ n and P{2). 

We will be using the following coordinatization of the Jacobi matrices 

00 00 
Q(l) = /+ + ££; a, Q(2) = I + + ^2^b l {i)E i ^ l (2.17) 

i 1=0 i 1=0 

One can immediately see that, for example, 



Q+(l) ).=S j ^ +1 + a (i)S iJ , (Q_(2) ),=R i S j ^i (2.18) 

where R4+1 = /tj+i//tj. As a consequence of this coordinatization, eq,( [2.14 ) gives in 
particular the important relation 

d 2 

———lnZ N (t,g) = Rn^ (2.19) 

5*1,10*2,1 

Finally let us quote from || the equation 

d 2 

InRj = R j+1 - 2Rj + Rj-i (2.20) 



5*1,1<9*2,1 

which justifies the name given to the hierarchy. 



2.1 H^i+oo constraints 

The W1+00 constraints (or simply VF-constraints) on the partition function for our 
two-matrix model were obtained in |nj,|J], by putting together both coupling con- 
ditions and consistency conditions. In other words the VF1+00 constraints contain all 
the available information. They take the form 

W^Z N (t,g)=0, W^Z N (t,g) = r > 0; n > -r, (2.21) 
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where 

Wj? = (-gYC^(l)-& + n n] {2) (2.22a) 
WW = (- 5 )"£M(2)-£ [ !r ] (l) (2.22b) 

The generators £^(1) are differential operators involving TV" and t\ while £^(2) 
have the same form with replaced by i2,fc- One of the remarkable aspects of ( [2.21]) 
is that the dependence on the coupling g is nicely factorized. The TVi+oo algebra 
satisfied by the £n^(l) has been written down in ref.Q,^. In general we have 

[41(1),4J(1)] = (sn - rm)£ [ :^- 1] (l) + . . . , (2.23) 

for r, s > 1; n > — r, m > — s. Here dots denote lower than r + s — 1 rank operators. 
The algebra of the £n^(2) is just a copy of the above one, and the algebra satisfied 
by the Wn and by the Wn is isomorphic to both. 

There is a sometimes simpler way to write the above generators. It consists in 
introducing the U(l) current 

OO tlj OO o 

J(z) = £ rt r z^ + - + E ^" V ( 2 - 24 ) 

r=l r=l r 

and defining the density 

r i 1 / \ n+1 

£N ( , ) = __ : (_ a + J(,)) : .l 
Then can be recovered as 

= Res z=0 {6 n \z)z n+k ) 

The above definition holds for both the 1 and 2 sector. One can also consider the 
fermionized version of the above formulas. This leads us to the Wi+oo constraints 



suggested in ref.[12| for the c = 1 string theory. One easily realizes, by using either the 
bosonic or the fermionic representation, that the latter are a subset of the constraints 
( |2.21 ). In fact they coincide with the cases n = —r and g = — 1 and can be written 



explicitly in the form 

d 



dti : r 

d 



Z N = £ [ : ] r (2)Z N (2.25) 
Z N = &\{l)Z N (2.26) 



dt 2 

This is our first link between two-matrix model and c = 1 string theory. 



3 Correlation functions of discrete tachyons 



In this section, from the W constraints (2.21) we calculate a subset of the correlations 
functions of two-matrix model in a very simple small phase space and identify them 
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with the correlation functions of the discrete tachyons of the c = 1 string theory. 
From now we set g = — 1. 

To start with let us introduce some simplified notations: 

h,k = tk, *2,fc = s k 

TrM? = r k , TtM* = a k 

The correlation functions of the two-matrix model are defined by 

d d d 9 

< T kl ■ ■ ■ Tk n (T h ...0\ m >= -K— ■ ■ ■ -K—^ • • • p mZ N (t, S) 

ot kl dt kn dsi, dsi 



Let us write some of the H^-constraints in this new language. W^\z^ = and 



W^Ii-^v = become, respectively 

oo 

Y kt k < T k-i » +Nti- < o-i >= (3.27a) 

k=2 

oo 

^ ks k < a fc _i > +iVsi- < n >= (3.27b) 

k=2 

Instead W^Zjv = and W$ ] Z N = give rise to the same equation 

oo oo 

kt k < r fc >= ^ ks k « (T k » (3.28) 
fc=l fc=l 

and so on. 

The TV constraints considered so far are exact, they contain contributions from 
all the genera. For simplicity we limit ourselves in this section to the genus contri- 
bution. The W constraints assume, in this simplified expression. To find it 
one can follow the homogeneity analysis of Q. Equivalently one can define x = N/(3, 
consider the rescalings 

p p p 

and keep the leading term in (3 in ( 2.21| ). Here /3 is an arbitrary positive number and 
is the genus part of lnZjy. Once this is done one quickly realizes that the genus 
version of ( p. 27a ), ( 3.27b| ) and (3.28) remain (formally) the same, except for the fact 
that N is replaced by x and now <C • 2> denotes only the genus contribution, not 
the complete correlation function. The latter simplifying convention will be followed 
until the end of this section. However the higher W constraints in general change 
form when reduced to genus 0. For example, W_ 2 Zn = becomes 

oo oo I— 3 

h,h=l 1=4 k=l 

oo 

+ 2x Y Hi < n-2 » +2x 2 t 2 + xt\ = <C a 2 » (3.29) 

1=3 
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Next we specify the small phase space in which we want to compute our correlation 
functions. This is the simplest possible one: t k = = s k , Vfc. We will denote the 
correlation functions calculated in such restricted parameter space as < • >, instead 
of < • >. 

Now we can set out to calculate the correlation functions. First of all let us 



differentiate both sides of ( |3.28| ) with respect to t kl , . . . , t kn and , . . . ,ti m . We get 

(kt + . . . + k n - h - ■ ■ ■ - l m ) < Tk x ■ ■ ■ Tk n a h . . . ai m >= 

which means that the correlation functions are nonvanishing only when 

fci + . . . + kn - h ~ ■ ■ ■ ~ lm = (3.30) 

We remark that this can be interpreted as a charge conservation. Next from ( |3.27a| , 
3.27b , |3,29|) we get immediately 

< a\T\ >= x, < 02T2 >= 2x 2 

In general from ( |2.25| , 2.26| ) we get 

<r h ai >= lx l 

Next one differentiate ( 3.27a ) w.r.t. t\ and si—%, and ( |3.29D w.r.t. ti,si—2, an d obtain 

< atai-xTi >= I < Ti-xai-x > 

< 0201-2TI >= 2x1 < (Tl-2Tl-2 > 

and so on. From this it is easy to conclude that 

< auo^kn >= kl(l - k)^- 1 (3.31) 

Proceeding the same way it is not difficult to arrive at 

< r k a h a l2 ...a lm >= khh • • . l m (k - l)(k - 2) . . . (k - m + 2)x k ~ m+l (3.32) 

where k = l\ + 1% + . . . + l m , m > 3. If we rescale t k — > kt k and s k — > ks k , the new 
correlation functions (3.32) become 

< r k a h a h . . . a lm >= (k - l)(fe - 2) . . . (k - m + 2)x k ~ m+l (3.33) 

These expressions are pretty familiar to those who are acquainted with the literature 
on the c = 1 string theory. If we set 

T k = r k , T_ k = a k (3.34) 

and interpret x as the cosmological constant, the correlation functions ( |3,33 ) and the 
like are nothing but the correlation functions of the discrete tachyons T k calculated 
in terms of the cosmological constant alone, [ [Df] , 

This is our second link between two-matrix model and c = 1 string theory. 

We remark that one can obtain many more results beside the genus correlation 
functions in the small phase space given above. One can explicitly calculate corre- 
lation functions in higher genus and in a larger space of parameters. However, since 
in this letter, our concern is to motivate the connection between two-matrix models 
and c = 1 string theory, we postpone these calculations to a future publication. 
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4 The dispersionless Toda lattice hierarchy 



We have already pointed out that the integrable hierarchy that appears in our two- 
matrix model is the Toda lattice hierarchy (see also Q). In two recent papers, 
|| and P], it has been suggested that the dispersionless Toda hierarchy underlies 
the Landau-Ginzburg formulation of the c = 1 string theory. To those familiar 
with integrable hierarchies it is almost evident that our integrable hierarchies ( |2.10b , 



2.10c , p7l2b| , 2.12c) coincide in the dispersionless limit with the dispersionless Toda 



lattice hierarchy of |2| and |3[. It is also apparent that the constraints introduced ad 
hoc in the latter references are nothing but our coupling conditions (|2.8| ), which are 
natural relations in the framework of two-matrix models. 

However, for most these are rather technical issues. Therefore we will spend this 
section to explicitly show the identification we just claimed. As our first step we 
rewrite the formulas of section 2 in the dispersionless limit. This can be achieved 
with a continuum limit. We showed in a series of papers (see, for example [H], ||), 
that there is actually no need to take this continuum limit, we can obtain the same 
results more neatly with the discrete approach, especially when all the genera are 
involved. However here we have to make a comparison with a continuum formulation. 
Therefore in this paper we will shift to it. 

The continuum dispersionless limit is obtained by promoting the matrix index n 
to a continuum variable x and by introducing the conjugate variable £ 

{(,x} = ( (4.35) 

We recall that the variables x and C, can be traced back to the discrete matrices 
p = J2n n E n ,n and I + || , respectively, and the Poisson bracket (|4.35| ) is nothing but 
the continuum version of the commutator 

[I+,p]=I+ 



In this limit the 'fields' ai(n) and 6/(n) on the lattice (see eq.( |2.17 )) are mapped 



into fields which are function of x (beside the coupling constants). Therefore we are 
going to have the following replacements 

oo 

Q(l)^L = C + X>rC r (4.36) 

r=0 

oo 

Q(2)^L = ( + Y / brC r (4.37) 



Similarly 



Moreover 



r=0 



P(l) -» M = rt r L r - x + XL' 1 + L ~ r ~ X (4-38) 

r=l r=l 

P(2) -^M = J2 rsM- 1 + XL' 1 + £ ^— Z^ 1 (4.39) 

r=l r=l ' 



oo oc i -\ 

lnVEV In* = J2 + ^ lnA i " E ^7 *T^ (0) ( 4 ' 4 °) 
r=l r=l 1 
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Now we have to introduce the continuous analog of the operation that maps a 
matrix Q to the matrix Q, which we defined at the beginning. This is the operation 
a which maps 

a(0 = | (4.41) 
and leave unaltered all the other quantities. Here R is the continuous limit of the 



'field' R n which we have defined after eq.Q2.18) 



The form the coupling constraints Q2.8| ) take in the dispersionless limit (remember 
that g = — 1), is 

M = a(L), L = a{M) (4.42) 

respectively. 

We can now write down the dispersionless versions of the hierarchies ( [2.10b , 2.10c , 



wr lL +' Lh £= {l ' m:) <443) 

dM dM 

_ = {i ;,M}, ^ = {M,ML} (4.44) 

Here the subscript + means the part of an expression contaning non-negative powers 
of C, while the - indicates the complementary part. Isomorphic hierarchies can be 
obtained by applying the a operator to these equations. 

It is now easy to make a comparison with refs.[|| and || and verify that the 
hierarchies are the same. To be more precise the correspondence of our paper with, 
for example, ref.|| is (our notations are on the left hand side) 

L «-> L, M <-> ML- 1 

M ^ ML, a{L)^L- 1 

Moreover the constraints (23), (24) of 0] are nothing but ( f4.42j ) above, and so on. It 
is apparent that the definitions and constraints introduced ad hoc in refs.||] and 0] 
are completely natural in the framework of the two-matrix model. 

This is the third element in common between the two-matrix model and the c = 1 
string theory. 

It has been suggested that the c = 1 string theory can be given a topological 
Landau-Ginzburg interpretation, ,|| , Q , . In our language the Landau-Ginzburg 
potential proposed by these authors coincides with M. In fact let us define 

, dM , dM dM 

0O= ^' - = ^- (4 - 45) 

We notice that if ij = = Sj Vz, j, then 

M = ~, L = C = Xi 
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and 



If we define 



= iC \ <p = <f>-j = jx j ( j 1 



^ j. j. j. ^ d 0i0i^fe 



where M' = Q^-, we find immediately that these three point function coincides with 
( |3.31 ). After a trivial redefinition the </>j's are the Landau-Ginzburg representatives 



of the % if M is assumed to be the corresponding potential. 



5 Discussion 

One could add other elements to the correspondence between the two-matrix model 
and the c = 1 string theory. A short example is the following. Eq. fl2.20 ), in the 
dispersionless limit, reads 

r) f) 2 7? 

ln#=^ (5.46) 



dt\ds\ d 2 x 
On the other hand, from eq.( |2.19| ) we have 

F) 2 

R = 6 J*® 
dt\dsi 

But the RHS of this equations is nothing but the correlation function calculated in 
the previous section, namely < T\C\ >= x. Putting all this together we obtain 

8 2 

u (0) = t^t^ (0) = lni? = hvx (5.47) 
ox 2 

This is the expected behaviour of the chemical potential in terms of the cosmological 
constant in the c = 1 string theory in genus 0. 

Even though we have not discussed the subject thoroughly (in particular for lack 
of space we postpone a discussion of the discrete states), we think we have given 
enough evidence that a two-matrix model underlies the c = 1 string theory. This 
claim may sound a priori surprising, since we are familiar with the idea that the c = 1 
string theory is described by a model of one time-dependent matrix. However it is 
more plausible than it seems at first sight, if one thinks that the connection between 
two-matrix model and c = 1 string theory is at the level of topological degrees of 
freedom. From this point of view let us recall that in [|J we showed two matrix- 
models contain all n-KdV hierarchies and the relevant A n series of topological field 
theories. It has been suggested recently fl4| that they might contain also the D series 
of topological field theories and, perhaps, other series. We do think that the two- 
matrix models are a perfect framework for a large family of topological field theories 
and will describe the topological field theory content of it in a forthcoming work. 
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